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Abstract 



When one tries to simulate quantum spin systems by the Monte Carlo method, 
often the 'minus-sign problem' is encountered. In such a case, an application of 
probabilistic methods is not possible. In this paper the method has been proposed 
how to avoid the minus sign problem for certain class of frustrated Heisenberg mod- 
els. The systems where this method is applicable are, for instance, the pyrochlore 
lattice and the Ji — J2 Heisenberg model. The method works in singlet sector. It 
relies on expression of wave functions in dimer (pseudo) basis and writing down the 
Hamiltonian as a sum over plaquettes. In such a formulation, matrix elements of 
the exponent of Hamiltonian are positive. 



The systems considered in this paper are Heisenberg models for spin one-half, with 
competing antiferromagnetic (AF) interactions - i.e. frustrated ones: 



where s is a vector of Pauli matrices, A C Z*^, Jij > (AF case). 

Frustrated systems are very interesting and hard to analyse and understand, 
both in classical version and especially in the quantum case. The source of these 
difficulties traces back to the large ground-state degeneracy in the classical version. 
The prototype of such system is AF Ising model on triangular lattice [T]. Such 
systems are very sensitive to perturbations. A consequence is a possibility of very 
complicated ground-state and finite temperature phase diagrams. This situation 
takes place, for instance, in an ANNNI model [2] (infinite number of phases, devil's 
staircase, etc). Besides of numerous efforts and important results [2], [3], [4] (for 
reviews, see [5], [6], [lO]), full treatment of such systems is not worked out so far. 

The situation for quantum frustrated antiferromagnet is even less clear. It is 
generally suspected that - in the case of strongly frustrated systems - the ground 
state emerging as a linear combination of many classical configurations is a feature- 
less, "spin liquid" state, i.e. state without long-range ordering, where correlation 
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functions fall-off exponentially [8]. However, one cannot exclude another scenario: 
"order by disorder" - exotic orderings absent in a classical version of these models. 
Such scenarios are moreover sensitive to the underlying lattice structure. (For a 
review, see for instance [9]). To my best knowledge, no general definite conclusions 
have been obtained so far. 

Among frustrated lattices, perhaps the most popular ones are: triangular; kagome; 
pyrochlore; square lattice with 'crossing bonds' (called also the Ji — J2 model). This 
last case is particularly interesting due to its possible relation with high-temperature 
superconductivity: Quite often one considers the t — t' two-dimensional Hubbard 
model as a 'minimal model' for HTSC [7]. Behaviour of this last model is still not 
fully understood. A natural starting point in such a study is the limiting case: half- 
filling and large coupling constant; under these conditions, the t — t' Hubbard model 
simplifies to Ji — J2 Heisenberg model. 

As a sample of natural questions in the study of frustrated systems one can men- 
tion the following ones (for definiteness, let us concentrate on the Ji — J2 Heisenberg 
model): 

• Nature of ground state: For which range of values of the ratio a = J2/J1 
we have an antiferromagnetic (Neel) ordering? Is there a spin-liquid state for 
strong frustration? 

• Describe the nature of crossover between ordered and disordered state upon 
increasing frustration. 

Exact results on the area of frustrated models are rather rare. For particular 
forms of interactions, there exist exact results for ground states, obtained by AKLT 
[TT] as well as related results [12]. It is however unclear if they can be generalized 
to more general forms of interactions. Some general properties of frustrated systems 
have been obtained in [14], [15], [I^ (they are important in the context of this 
paper). One can mention also quite a few approximate reliable results, for instance 
[T3] (based on BCS-like ansatz on wave function). 

One of general tools used to calculate the partition function Za and thermody- 
namic functions for quantum spin systems is an application of Lie-Trotter product 
formula. Let us describe the general setup of certain version thereof, i.e. the Suzuki 
approach [T7] . 

Lie-Trotter product formula [TT], [IH] states that if A, 5 - finite-dimensional 
matrices, then 




Using this formula, one calculates Za in the following way: 



1. Write: 



H = Hi + H2, 



in such a way that Hi, H2 are sums of commuting operators. 
2. Using Lie-Trotter formula, we have: 
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■ {0i2n-l\e " |a2n)(a2n|e " |ai) (1) 
(here {aj} is a basis in the Hilbert space of system states). 

If it happens that above matrix elements are positive, then the hfe is easier, as 
one can apply probabilistic techniques: 

• Monte Carlo method - in numerical aspects; it is called the quantum Monte 
Carlo [20], 

• contour expansion techniques [18] or "stochastic geometry" methodology in rig- 
orous studies; as examples, one can mention spin chains [21], or Bose-Hubbard 
models [22]. [23], 

For numerous important cases, matrix elements are positive. It is the case, for 
instance, of the quantum Ising model in transverse magnetic field, ferromagnetic 
Heisenberg model, XY model, Falicov-Kimball model (for a review, see [20]). These 
systems as well as numerous other ones have been successfully studied with the use 
of quantum Monte Carlo method. 

But alas! In general, matrix elements are not positive. (This is famous "minus- 
sign problem" in the quantum Monte Carlo method). 

In certain cases, this problem can be overcomed. For instance, if one considers 
antiferromagnetic Heisenberg model, then for simplest choice of basis states in ([T]) 
(Ising basis), the problem is present. But it can be overcomed for the model on 
bipartite lattices, using more sophisticated techniques [2l]. This is also the case of 
the Hubbard model on bipartite lattices and for half-filling [25]. Some results for 
frustrated antiferromagnets have been reported [26]. However - to my best knowl- 
edge - the solution of the minus-sign problem is still lacking for general frustrated 
antiferromagnetic Heisenberg models. 

The goal of this paper is elaboration of the quantum Monte Carlo scheme for 
certain class of frustrated Heisenberg models. Using this scheme, matrix elements 
obtained are positive. 

This scheme concerns the Ji — J2 Heisenberg model and holds under certain 
conditions. There are: 

1. Presence of reflection symmetry in the system, 

2. We restrict ourselves to the singlet sector of the system (i.e. we assume that 
the total spin of the system is zero). 

The first assumption is not too restrictive and is technical only; it is necessary to 
apply the Lieb-Schupp theorem (discussed below). The second one is more serious. 
But we can argue as follows. For certain class of frustrated antiferromagnetic Heisen- 
berg models (including the pyrochlore and Ji — J2 model but not the triangular or 
kagome lattices), we have Lieb and Schupp theorem [ll],[T^, [16] stating that the 
ground state of such systems is singlet. One then can hope that performing the MC 
simulation in the singlet sector at finite temperature T, and then tending with T to 
zero, we will obtain reliable properties of ground state of such a system. 

The sketch of the scheme is as follows. It is well known that positivity of matrix 
elements is a problem of the choice of basis in the set of 'intermediate states' |aj)(aj| 
in formula ^ (instructive examples can be found in [2^). If one chooses the basis 
being a tensor product of Ising states (it is perhaps the simplest choice) as the basis 
of 'intermediate states', then some of matrix elements of Hamiltonian are negative. 

Assume however that we work in the singlet sector. It is known that every singlet 
can be built up from dimers, i.e. two-spin wave functions of total spin equal to zero. 
(This will be discussed in the Section [2l). Such a form of singlets have been used in 
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numerous papers: [8], [27], [28], [29]. Consider now the system defined on the 1? 
lattice or, more generally, on a bipartite one (i.e. composed of two sublattices called 
A and B in such a way that only neighbours of A type sites are B type sites and vice 
versa). Moreover, let us impose the condition that one spin of every dimer belongs 
to A sublattice and the second one to B sublattice. Now, consider the model with 
nn interactions (or, more generally, with unfrustrated ones). It turns out that (for 
details see below Sec. [3l) that matrix elements of the Hamiltonian are non-negative. 

But if we consider the Heisenberg model with frustrated interactions (for defi- 
niteness, take the Ji — J2 Heisenberg model), then some of matrix elements are still 
negative. How to cure this situation? The idea is as follows: Let us write the Hamil- 
tonian as a sum of the plaquette terms, i.e. four-spin Hamiltonians being defined 
on the 2x2 plaquettes on the lattice. One can hope that in such a situation, neg- 
ative contributions coming from nnn interactions will be compensated by positive 
ones coming from nn interactions. It turns out this is the case: Matrix elements 
calculated with the use of plaquettes are positive. The calculations are presented in 
Seed 

The paper is organized as follows. In the Sec. [21 the general setup is introduced: 
construction of singlet wave functions from dimers is explained, and the scalar prod- 
uct of two singlet wave functions is calculated and interpreted in geometrical terms. 
In Sec. [3l the matrix elements of the Hamiltonian written as a sum of two-spin 
interactions are calculated. 

The central part of the paper is SecHJ where matrix elements of the Hamiltonian 
written as a sum of four-spin (plaquette) interactions are calculated. Moreover, 
the (not quite obvious) positivity of matrix elements is proved there. The Sec. [5] 
contains summary and discussion of results obtained as well as considerations how 
to generalize the results. In the Appendix, technical tools used in calculations are 
presented (spectral resolution of self-adjoint operators and its application to two- 
and four-spin hamiltonians). 



2 Dimers, singlets and all that 

Consider the dimer, i.e. the singlet wave function localized on sites 

Assume now that the total number of spins is even, i.e. we are dealing with 2A'^ spin 
system. Then every singlet wave function ^'2Af can be built from dimers: 

^2Af = ^Ci^,...,iN;ji,-,jNihji){i2j2) ■ ■ ■ (Wiv), (2) 

[27] . [28] . [29) . This representation is non-unique for > 1 (the set of all such dimer 
products is an overcompleted set of vectors spanning all the space of singlets). 

Consider now the square lattice. It is bipartite one, and all considerations below 
refer also to such lattices. Divide the lattice into two kinds of sites: A- and i?-type 
sites. We demand that in jj]) 

ik G A, jk € B for all A; = 1, . . . , A^. (3) 

Also in this case the set of all dimer products is an overcompleted set (for N > 2) 
in the vector space of singlets. 

Consider now some singlet wave function on the lattice, which is a product of 
dimers. Such a function possess a natural geometric interpretation [2^. Every dimer 
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Figure 1: An example of dimer wave function on 4 x 4 lattice. Two sorts of sites of the 
bipartite lattice are represented by circles and heavy dots. Dimers are represented as 
'bonds' linking lattice sites of opposite kinds. 

i'ikjk) can be illustrated as a 'bond' linking lattice sites and jk (remember ik S A, 
jk G B). Notice that every lattice site is occupied by the end of exactly one bond; 
in the other words, dimers are 'closely packed'. Such a situation is illustrated on 
Fig.m 

Consider now the scalar product of two such functions ^'i, ^'2: 

Sij = {^i\^2) (4) 

Let us draw both functions on a lattice. Such a situation can be viewed as a 
set of closed polygons. Every such a polygon is formed by dimers belonging to 
^'i, ^'2, ^'i, ^'2; • • • (in alternating manner); the number of bonds forming this poly- 
gon is even. It is illustrated on the Fig. [2l 

Consider first the situation where wave functions ^'i,\l'2 correspond to single 
non-trivial polygon on the lattice. (We call the polygon nontrivial if it is not a 
'double bond', i.e. if its length is 2L, L > 1). It is a matter of straighforward 
calculation (it follows also from the 'reduction principle', see below) to show that 

m 

(^11^2) = (5) 

This result can be generalized to the situation where ^'i, 'I'2 correspond to family of 
polygons: Pi of length 2Li, . . . , of length 2Lfc. In such a case, we have |27| : 

(^i|^2) = 2-E-i(-^'-^)l (6) 



3 Matrix elements of exponens of the Heisen- 
berg operator: two-spin form 

In this Section, we consider the model with nearest-neighbour interactions. (Al- 
though the main interest of this paper are frustrated models, considerations of this 
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Figure 2: Two dimer functions dashed line and "^2' continuous line) and polygons 
formed by them. Every site is occupied by exactly one end of dimer belonging to and 
the same for '^2- On the picture there are two trivial polygons (formed by 2 dimers) and 
three non-trivial ones (formed by 4 dimers). 



section can be treated as a warm-up and a presentation of techniques used in the 
next section). The Hamiltonian is of the form 

H = jY,hij, (7) 

(ij) 

where hij = Sj • sj. 

Write an expression jT]) in the form 

H = Hi + H2 + H-3 + H4, (8) 

in such a way that Hi, . . . , H4 are sums of commuting operators. One possible way to 
achieve this goal is as follows. Every site index i is in fact a two-index: i = {ix,iy), 
where ix is horizontal index on the lattice and iy is vertical index. Divide the 
Hamiltonian ^ into Hi, . . . , H^ in the following way: 

• i^i is a sum of these operators hij where i,j are of the form: {2k, 21), {2k, 2^+1). 
Denote it as type 1; 

• -ff2 is a sum of /ij^'s where i,j are of the form: {2k, 21), {2k + 1, 21) - type 2; 

• i^s is a sum of h'-s where i,j are of the form: {2k, 21 + 1), {2k + 1, 2/ -|- 1) - 
type 3; 

• Hi is a sum of /ij^'s where i,j are of the form: {2k + I, 21), {2k + 1, 2/ -|- 1) - 
type 4. 

Consider now a matrix element (^/l exp(i('iffc)|^j), k = 1,...,4 in order to 
check its positivity. Every operator H^ is a sum of commuting operators, so, if 
Hk = hij, then ex.p{KHk) = ex.-p{Khij) - so one can write: 

i,j of type k i,j of type k 

{^i\eMKHk)\^j) = 
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(here M is a number of operators hij in i^fe). We can make a conclusion that if 
matrix elements of the operator exp{Khij) is positive, then the matrix element of 
exp{KHi.) is also positive. 

Let us calculate the matrix element of two-spin operator hij: 

{^i\eMKhj)\^2) (9) 

where |^'i),|^2) are dimer functions, K = —(3J/N. Notice that antiferromagnetic 
case J > imply K g] — oo, 0[. 

Below, the following notation will be useful: 

= exp(i^/2), £3 = exp(-3X/2). 

Notice that G]0, 1[, es g]1,oo[. 

Let us consider first the situation where i,j are nearest neighbours. 
We have three sorts of situation: 

a) The operator exp{Khij) acts on which contains the (ij) dimer; it is 
illustrated on Fig. [SK- The value of the matrix element is 

(^-21 exp{Kh,,)\^i) = e3(*2|^'i>; (10) 

we see that this element is positive. 

b) The operator ex.p{Khij) is localized on a bond connecting two different poly- 
gons; see Fig. [Sb- We have: 

(^2| exp{Khij)\^i) = i(3ei + e3)(^'2|'f i) (11) 

which is also positive. 

c) The operator exp{Khij) acts inside one connected polygon, but there are at 
least three polygon 'bonds' between i and j sites; see Fig. [3b. In such a situation we 
have: 

(^-21 exp(K/i,,)l^i) = e3(^'2|^'i); (12) 

clearly, it is also positive. 

Similar matrix elements (but of the operator hij instead of its exponent) have 
been calculated in |27j. In the case considered here, they have been derived by 
straightforward calculation with the use of simple algebraic tools: spectral resolution 
of the operator hij and the 'reduction principle'. Details of calculation of matrix 
elements is presented below. 

A matrix element for situation a) (presented on Fig. [3h) can be calculated im- 
mediately. Namely, in this situation the operator ex.p{Khij) acts for dimer function 
(ij), which is an eigenfunction of hij (and of course of exp{Khij)). Using property 
([24|l we get the expression (fTO]l . 

For situations b and c) (presented on Figs. [3b and [3]:), we apply the reduction 
principle first. Such a reduction is performed in two steps R.I and R.II. 

R.I We examine the quotient q: 

_ (^2|exp(K/i,,)|^i) 

^'i and ^'2 functions form closed (not necessarily connected) polygon P. Assume 
that these functions can be written as: ^'i = ^i{lk){nm), ^'2 = ^2il'n^)inp). Let us 
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Figure 3: a), b), c) - three kinds of matrix elements of the operator exp{Khij) listed 
above. The action of this operator is denoted as a bold line; gray and dashed lines denote 
dimers forming the functions and \l/2, respectively, c') is obtained from c) by one 
step of reduction principle R.I; upon this operation, neighbouring dimers Di and D2 are 
eliminated, b') is obtained from b) by three steps of reduction R.I. c") has been obtained 
from c') by two steps of R.II. 
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remove the site / and identify k and m sites (this situation can be viewed as removing 
dimers Di = (kl) from and D2 = {mk) from ^'2 with subsequent identifications 
of sites k and m). This way, we obtain 'reduced' functions ^'f = <I>i(n/c) and 
3'^ = ^2{np) and 'reduced' polygon P^. We assume that all sites k,l,m,n,p are at 
a distance greater than 1 from i and j. Then, the quotient q obtained for functions 
and ^'2 ([13]) is equal to the quotient calculated for functions and 

o (^f|exp(irfe,,)|^f) 

Less formally, we can say that the quotient q will not change if we remove two 
neighbouring dimers from the polygon P. 

Proof. Let us remember how to calculate the scalar product for dimer wave 
functions: After expressing them in the 'plus-minus' basis, we sum over all sites and 
spin degrees of freedom for every site. 

Let us consider the matrix element (^'2] exp(i('/ijj)|^i). Let us expand dimer 
functions (Ik) , (nm) , (Im) , (np) in the base of 'plus-minus' functions. Then, let us 
sum over spin indices of sites I and m. After straightforward calculations, remem- 
bering about normalization factors for dimers and using orthonormality relation of 
spin functions on arbitrary site r: (ro-|ro-') = (^o-o-', we obtain: 

(M/2I exp{Khij)\^,) = I exp{Khi,)\^f) (15) 

Analogous calculation gives 

(*2|^'i) = li^^l^f) (16) 

so we obtain desired equality q = qn- 

After (possibly, multiple) application of reduction principle to situations as in 
Figs.Eb andEb, we get situations such as on Figs.Eb' and Eh'. This way, reduction 
principle makes possible the calculation of matrix elements for situations, where 
wave functions occupy at most 6 sites. 

R.II By straightforward calculation one obtains that both cases illustrated on 
Figs.[3l:' and[3h" give equal value of q. (This step is not necessary, but its analogon 
for plaquettes will be useful due to oeconomy reasons). 

Two cases obtained after reduction, i.e. the ones illustrated on Figs. [3b' and[3l3' 
(or, equivalently, [3):") can be calculated immediately. The line of calculations is as 
follows: For the wave function ^'i, one passes from dimer form to the 'plus-minus' 
basis. Then, one calculates an action of the operator exp{Khij) on making 
use of spectral resolution of the operator hij given by ([25]) . ([26ll as well as formula 
([^ . And last, one calculates the scalar product of expression obtained with the ^'2 
function, expressed in the 'plus- minus' basis. 

We can conclude this section by an assertion that matrix elements for nearest 
neighbour Heisenberg model are positive, so we have no minus sign problem here. 
In my opinion, such a result can be viewed as an interesting one, but not exciting: 
There are other approaches, where minus-sign problem has been overcomed [24] . 

m, m- 

Now, let us consider the frustrated case (i.e. the Ji — J2 Heisenberg model). In 
this case, however, the matrix elements are in general not positive. As an example, let 
us mention situation analogous to the case b) above, but where i,j are next-nearest 
neighbours. In such a case, we have 

(^2| exp{Khij)\^i) = -^(3e^ + e3)(^'2|^'i), (17) 
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1 2 




Figure 4: Elementary plaquette on which the plaquette hamiltonian (fTSl) is defined. Sites 
are arranged anticlockwise. 



which is not positive. In the other words: For frustrated ( Ji — J2) model, where both 
nn and nnn interaction are present, the minus-sign problem still exists. 

How to cure the problem? The idea is as follows: Write the Hamiltonian as a 
sum over plaquettes (4-spin) sets. One can hope that negative contribution will be 
compensated by positive one. (It is not obvious a priori, as two-body operators 
entering into plaquette operator does not commute in general). It turns out that 
in such a formulation, the matrix elements are positive; details are described in the 
following Section. 

4 Matrix elements of exponens of the Heisen- 
berg operator: plaquette form 

Consider the Heisenberg model on a (subset of) square lattice. We assume that 
there are nn and nnn interaction. For concreteness, we consider the Ji — J2 model, 
but all considerations apply also in the case of 'pyrochlore' lattice and some others. 
We assume that the system exhibits the reflection symmetry (remember in such a 
case, the ground state is singlet). 

The setup for wave functions is the same as previously: We assume that wave 
functions are built up from 'bipartite' dimers. 

Consider the Hamiltonian defined on a square plaquette: 

hn = hi2 + /l23 + ^34 + h4i + a{hi3 + /l24), (18) 

(see FigHI), where: hij = Sj • Sj, a = J2/J1. 
We have: 

H = Y,hD,^ (19) 

igA 

Analogously as before, let us write the Hamiltonian in the form 

H = Hi + H2 + H3 + H^, (20) 

where each of terms Hi, . . . ,Hi is a sum of commuting plaquette operators. This 
division can be done, for instance, in the following form. Plaquette index in (fT^ is 
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Figure 5: Illustration of the second step reduction (R.II.plaq) action, a) corresponds to 
initial configuration; b) is the configuration after single action of (R.II.plaq); on c), after 
three further actions of (R.II.plaq). All three configurations have the same value of q\j. 



in fact a two-index: 



H 



• ,iy), where ix - horizontal component, iy - vertical one. 

, Hi are defined as: 

• i^i is a sum of these operators /in^j where i is of the form: {2k, 21). Denote it 
as type 1; 

• is a sum of these operators /in^j where i is of the form: i2k + 1, 21). Denote 
it as type 2; 

• i^a is a sum of these operators /i^j where i is of the form: {2k, 2Z + 1). Denote 
it as type 3; 

• Hi is a sum of these operators /i^ j where i is of the form: {2k + 1,2/ + 1). 
Denote it as type 4. 

As in the previous Section one shows that if matrix elements of elementary pla- 
quette operator (^'/| exp(i^/iG)|'I' j) are positive, then matrix elements (^'/l exp{KHk) 
are also positive. 

Let us calculate the matrix element of h\j- This will be done in two stages. In 
the first step we apply the 'reduction principle'. Its first stage R.I.plaq is identical 
as in the case of two-site hamiltonian R.I, i.e. it allows to decrease the number of 
neighbouring polygon edges by two without change of the quotient q\j, where q\j is 

(*2|exp(E:/in)|1'i) 



(^'2|^l) 



(21) 



The proof of R.I.plaq is also identical as previously, so we will not repeat it. This 
way, all calculations are reduced to cases where wave functions occupy at most 12 
sites. It is possible further reduction. This second step of reduction is similar to 
R.II but not identical: 

R.II.plaq: Assume that in the configuration there appear a square consisting of 
the following sides: One edge belong to the hamiltonian plaquette (say, this is {i,j) 
side); two sides are dimers belonging to the function \I'i: {i,l), {k,j); and the last 
side is dimer belonging to the function "^2- (let it be {k, I)). In such a configuration, 
one can eliminate two sides {k,j) and {k,l) and replace the square {i,j,k,l) by one 
bond {i,j) with dimer {i,j). 

The proof of R.II.plaq can be obtained by a straightforward calculation. An 
example of its action is illustrated on Fig. [5l 

This way, after reductions, we have to calculate matrix elements of the plaquette 
hamiltonian with wave functions defined on at most 8 sites. 
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Figure 6: Reduced configurations for plaquette Hamiltonian. 



It turns out that there are eight such configurations. They are illustrated on 
FigEl The matrix elements corresponding to them can be calculated in a straight- 
forward manner - similarly as in the previous Section for the case of two-spin Hamil- 
tonian. Basic steps of such calculations are: 

1. We express dimer functions in the 'plus-minus' basis. 

2. We calculate the action of the exp(K/in) operator on the function using 
the spectral resolution of the operator h\j (necessary formulas are collected in 
Sec. 16. 2|) together with the formula ([2^ . 

3. And last, we calculate the scalar product of the expression obtained above with 
the wave function ^'2 (expressed in the 'plus-minus' basis). 

Calculations are straightforward but lengthy, and they have been performed with 
the use of symbolic calculations programme (Maple). The results are summarized 
in the Table below. 

In formulas below, the following notation has been used y = exp(— K); a = a/2. 
We assume that < q < 1 (as for these values of a we have warranty that the 
ground state of the system is singlet [E], [31]). Remember that for < T < 00, we 
have 1 < y < 00. 
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Table: Values of matrix elements of the quotient for eight reduced configura- 
tions on Fig. 



No 




I 


T 7 TTT n 7. 7 T 7 T T"^ 

i iy6a ^ y2-a ^ g^a ^ g^i-a ^ ^y-l-a^^ 

— ^ T 7 TZ ^ Z Z T :T^\ 


II 


T 7^ !T Z T Z T r"^ 


III 




IV 




V 




VI 




VII 




VIII 





Some of expressions in the Table are evidently positive (there are: I, II, V). It 
turns out that remaining expressions are also positive; proofs are presented below. 
Remember that a G [0,^] and y g]1,oo[. Factors before expressions for are 
skipped. 

• Situation III. We have: 

_yZa ^ yl-a ^ ^yl-a ^ ^y-l-a > _y3a ^ y2-a ^ y^a^_-^ ^ y2-Aay^ 

but for 2 - 4a > and for y > 1 we have (-1 + y^""^") > 0. 

• Situation IV. 

and this is the same expression as in III. 

• Situation VI. 

SO we again obtain case studied in III. 

• Situation VII. 

y^" + 3y2-'^ - Qy^ + Qy^-^ > 6(-y« + y^-") = 6y"(-l + y^-^") 

and, analogously as in III, conditions: 1— 2a > 0, y > 1 imply that — l+y^"^"^ > 
0. 

• Situation VIII. This is the most complicated one. Write first: 

y3a ^ y2-a _ g^a ^ ^yl-a ^ ^y-l-a ^ y3a ^ y2-a _ ^ya ^ 5y-l-« + 3(-y'^ + y^-"). 

Due to argumentation identical as in III, the last term (in parentheses) is pos- 
itive: — y° + = y"(— 1 + y^~^") > 0. So, it is sufficient to show positivity 

of y^" + y2-a _ 3ya ^ 5y-l-a ^ ya(y2a ^ y2-2a _ 3 _^ g^-l_2a^)^ 

Consider two extreme cases (i.e. a = and a = ^) of the expression in 
parentheses above. We have for a = 0: 

l + y2_3 + 5y-i >y-2 + i + - = f^-^) + - > 0, 

y y \ VvJ y 

and for a = I 

y + y - 3 + 5y-2 = y-2(2y3 - 3y^ + 5) = y-\y + l)(2y2 - 5y + 5), 
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which is greater than zero for y > 1, as it fohows from elementary considera- 
tions. 

For intermediate values of a, i.e. < a < i, let us write 

y2a ^ y2-2a _ 3 ^ 5y-l-2a ^ y-l-2a^yl+4a ^ ^3 _ ^yl+2a ^ 5) 

and notice that > _yi+2a > — y^, so we have 

yl+ia ^y3_ 3yl+2a + 5 > y + ^3 _ 3^2 ^ 5 ^ ^ ^) (^2 _ 4^ ^ 5) > 

thus establishing positivity of VIII. 

Let us summarize this section by the statement that in the dimer basis, and for 
the Hamiltonian written as a sum of plaquettes, matrix elements of the exp{Kh\zi) 
are positive. 

5 Summary, conclusions 

The technical tool to study (certain class of) frustrated systems has been developed, 
which (hopefully) would allow application of probabilistic techniques. 

This paper was devoted to the Ji — J2 Heisenberg model in any dimension (and 
models whose Hamiltonians can are sums of frustrated plaquette hamiltonians; the 
pyrochlore lattice is perhaps the most typical example). It would be tempting to 
generalize the method to other frustrated systems. Such a generalization probably 
can be realized in the case of other systems exhibiting reflection symmetry, for 
instance 3d Heisenberg models with frustrated cubes; of course, one have to calculate 
the matrix elements of Hamiltonians for frustrated units. 

Generalization for frustrated systems exhibiting no reflection symmetry (such as 
kagome or triangular ones in two dimensions) is less obvious. For such systems, 
matrix elements of frustrated units can be positive (author calculated such elements 
for triangular lattice and positivity holds also in this case). But, on the other hand, 
the method relies heavily on the assumption that we are working in the singlet 
sector. It has been proved that the ground state(s) of the Ji — J2 Heisenberg model 
is singlet [Il]-[l6], [M], but for triangular or kagome lattices it is not known. (The 
answer is probably positive, but the proof - as far I know - is lacking). If the ground 
state is singlet, then also in the case of triangular lattice one can try to simulate the 
ground state using this method. 

The next point is the technical one: Full Monte Carlo simulation of dimer- 
like models is rather difficult task [32]. It is not clear how difficult would be an 
implementation of the actual method; this paper is devoted only to development of 
the scheme. However, certain attempts towards concrete computational realization 
of this method are in progress. 

It would be interesting to try to develop analogous method in the sector of S 
arbitrary, not only 5" = 0. If successful, it would be possible to perform qMC simula- 
tions in arbitrary temperatures, not only in low-T region (as in the present version). 
Moreover, it would be also possible to simulate systems where is no warranty that 
the ground state is singlet (for instance, the triangular lattice). There are some 
indications that the procedure described in this paper could be generalized for S 
arbitrary. However, at this moment it is too early to say something definite. 
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6 Appendix: Techniques used in calculation of 
matrix elements 



Let A be the self-adjoint operator in finite- dimensional Hilbert space TC, dimH = N: 
let Sp(A) = {Aj} - spectrum of A, Vi - the subspace corresponding to the eigenvalue 
Aj, dimVi = rii. We have: X^jWi = N,H = ®Vi . Every such operator A can be 

i 

represented in the form of the spectral resolution: 

A = Y,\Pi, (22) 

i 

where - orthogonal projection onto corresponding subspace Vi. These projections 
possess well known properties: YliP'i — {Idl-i ^ the identity operator in TL); 
Pf = Pi for every i; PiPj = PjPi = for i / j. Every such projection Pj onto 
eigensubspace Vi can be calculated from the famous formula 

m 

Pi = J2MM (23) 
k=l 

where is fc-th vector of orthonormal base spannning the Vi subspace. 
If f{x) - analytic function, then we have: 

/(^) = ^/(A,)P, (24) 

i 

For the sake of completeness, we present below expressions for projections, which 
appear in the spectral resolutions of operators hij and h\j. 



6.1 Spectral resolution of the operator hij 

For two spins, the state space Ti.2 is four-dimensional. The Hamiltonian hij = 
Si ■ Sj commutes with the total spin operator, so eigenvalues of hij can be classified 
accordingly to angular momentum quantum numbers. Eigenvalues of hij are: Eq = 
— I (total spin S is zero, i.e. the state is singlet; its multiplicity is one), and Ei = j 
(here total spin S is one, i.e. it is a triplet. We have three triplet states with z-ih. 
component of angular momentum equal to -M,0, — 1; all of them have the same 
energy, so the multiplicity of Ei is 3). 

Let us choose in 7^2 the standard ('plus-minus') basis ei, . . . ,64: 



ei 



\i+j+),e.2 = \i+3-),ez = \i-j+),e4 = 



In this basis, the Hamiltonian hij is given by the matrix: 



Si ' Sj 



10 
0-120 
2-10 
1 



(25) 



Projectors Pq (onto singlet subspace) and Pi (onto triplet subspace) are: 



Pn = 













1 




-1 

1 




Pi = 



2 

110 

110 

2 



(26) 
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6.2 Spectral resolution of the operator h\j 

The space of states for the system of four spins Ti^ is 2^ = 16-dimensional. In such 
a situation, it is again very useful to exploit properties of the angular momentum 
operator and classify states accordingly to spin quantum numbers. 

The plaquette Hamiltonian is given by (fT8]l . Remember a = J2/ Ji; it is more 
convenient to introduce a = a/2. The Hamiltonian (fTSl) commutes with total angular 
momentum operator S. For eigenstates and eigenvalues of the Hamiltonian, good 
quantum numbers are S (total spin) and M [z-th component of spin). Let us classify 
states accordingly to the value of M first, and then, in every sector with given M, 
classify the states accordingly to S. For each such state, we will give their energies 
and corresponding projectors. Classification of states with respect to values of M 
proceeds as follows. 

Among all 16 states, we have: 

• One state with M = 2 (plus twin state M = —2). Both of them share S = 2. 

• Four states with M = 1. Among them, there are: One state 5 = 2, M = 1 and 
three states 5 = 1, M = 1. There are also twin states for M = — 1. 



six M = states. There are: one state S 
0; two states S" = 0, M = 0. 



2, M = 0; three states S = 1,M 



We can consider only states with non-negative values of M, as all of them possess 
their twins for — M. 

Consider now all M sectors. Eigenvalues will be denoted as : e''^\j [n - index 
of the state), and corresponding projections by Pg^lf- 

• M = 2 sector. 



(2) 



This is subspace spanned by one base vector e\ 
in this sector is simply a number 1 + a, which is of course also the eigenvalue 



The Hamiltonian 



£"2^2- The projection is trivial. 
M = 1 sector. 

This subspace is four-dimensional. As a basis, let us choose: 



4'^ 



,(1) 



) ^4 



(27) 



The Hamiltonian in this basis is given by: 






1 


2a 


1 


1 





1 


2a 


2a 


1 





1 


1 


2a 


1 






(28) 



Eigenvalues (i.e. energies) and corresponding projectors are: 



(1) 

'2,1 



1 + a. 



PI 



(1) 



2,1 



1111 
1111 
1111 
1111 



(this is M = 1 component of quintet). The remaining three states are triplets 
(more precisely, the M = 1 components thereof). Two of them are degenerate: 
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the subspace spanned by eigenvectors of e[^1 is two-dimensional: 



E 



(1) 
1,1 



-a P(') 
a, r-^^^ 





1 





-1 





1 





1 





-1 


2 


-1 





1 










-1 





1 



E' 



(2) 



(2) 

1,1 



And last, consider 

The M = sector. The basis is: 



=(0) - 



-1 



+ + 



.(0) 



+ /'^3 



(0) 



1 


-1 


1 


-1 


-1 


1 


-1 


1 


1 


-1 


1 


-1 


-1 


1 


-1 


1 




+ - 








+ - 




^4 



.(0) _ 



=(0) _ 



+ 



The Hamiltonian in M = sector is: 



h = 



Eigenvalues and projectors are: 

-) The M = component of quintet: 



" -2a 





2a 


2a 


1 


1 





-2a 


2a 


2a 


1 


1 


2a 


2a 


-2a 





1 


1 


2a 


2a 





-2a 


1 


1 


1 


1 


1 


1 


-2 + 2a 





1 


1 


1 


1 





-2 + 2a 



Ei]l = l + a, 



2,0 



1 1 1 
1 1 1 
1 1 1 
1 
1 
1 



-) The M = components of triplets: 



-^1,0 ~ -'^1,0 





" 1 


-1 











" 




-1 


1 














1 








1 


-1 








2 








-1 


1 

















































(remember that subspace spanned by eigenvectors corresponding to E 
two-dimensional), and 



-^1,0 — 



-1 + a, 





" 














" 






















1 




















2 


































1 


-1 
















-1 


1 
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-) And last, singlets: 



-^0,0 



-3a, Pr 



(1) 



0,0 



1 1-1-10 

1 1-1-10 

-1-11 10 

-1-11 10 







(34) 



E, 



(2) 
0,0 



-2 + a, Po^5 = - 



1 

12 



1 1 1 1-2-2 

1 1 1 1-2-2 

1 1 1 1-2-2 

1 1 1 1-2-2 

-2 -2 -2-2 4 4 

-2 -2 -2-2 4 4 



(35) 
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